A spin-boson thermal rectifier 
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Rectification of fieat transfer in nanodevices can be realized by combining the system inherent 
anharmonicity with structural asymmetry, we analyze this phenomenon within the simplest anhar- 
monic system -a spin-boson nanojunction model. We consider two variants of the model that yield, 
for the first time, analytical solutions: a linear separable model in which the heat reservoirs con- 
tribute additively, and a non-separable model suitable for a stronger system-bath interaction. Both 
models show asymmetric (rectifying) heat conduction when the couplings to the heat reservoirs are 
different. 
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The heat conduction properties of nanojunctions at- 
tract attention for two reasons. First, heating in 
nanoconductors, a crucial issue for their operation and 
stability, is determined by both heat release and con- 
duction in such systems. Secondly, as with electronic 
conduction, the restrictive geometry raises fundamental 
questions concerning the relationship between transport 
processes in microscopic systems and their macroscopic 
counterparts. Indeed thermal transport properties of 
nanowires can be very different from the corresponding 
bulk properties as is demonstrated by the recent confir- 
mation p] of the prediction that the low temperature 
ballistic phonon conductance in a 1 dimensional quan- 
tum wire is characterized by a universal quantum unit. 
Also of considerable interest are studies that confront the 
macroscopic Fourier law, J = —KVT , a linear relation- 
ship between the heat current J and the temperature gra- 
dient VT that defines the thermal conductance with 
heat transport on the microscopic scale 0, S IE IS 13 ■ 
Harmonic chains were repeatedly discussed theoretically 
in these contexts and considerable experimental progress 
was also made For reviews see Refs. [lllfl^ . 

An intriguing mode of behavior often addressed in the 
study of transport devices is current rectification, allow- 
\ ing larger conduction in one direction than in the oppo- 
site one when driven far enough from equilibrium. Such 
' phenomena were extensively studied for electronic con- 
duction in molecular junctions, but much less so for ther- 
] mal nano-conductors. For a harmonic thermal conductor 
connecting (by linear coupling) two (left (L), right (R)) 
harmonic thermal reservoirs that are maintained at equi- 
librium with the temperatures and Tr, respectively, 
heat transfer is a ballistic process and the heat current J 
can be recast into a Landauer type expression 0, 0, 0| 

Jiuj)^ T{uj)[nL{uj) - nR{uj)]ujdLU, (1) 
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where T^uj) is the transmission coefficient for phonons of 

frequency cu and riKiuj) = (e^^'^ - l) ; = (kBT)^^; 
K — L,R {h = 1) are Bose-Einstein distribution func- 
tions characterizing the reservoirs. Obviously, this ex- 
pression is symmetric to interchanging the reservoirs tem- 
peratures and cannot show rectifying behavior irrespec- 
tive of any asymmetry in the system structure. 

In contrast, Terraneo et al l^j have shown numerically 
that rectifying behavior is obtained by replacing the inte- 
rior part of a classical harmonic chain by an anharmonic 
segment. An example of a similar behavior by a some- 
what simpler model is shown in Fig. ^ The model is 
defined by the TV-particle Hamiltonian 

N N-1 , 

1=1 i=l 

D (^e-"(^i-'') +D (^e-"(^-^") - l)^(2) 

supplemented by damping and noise terms operat- 
ing on particle 1 and N to simulate the effect 
of two thermal baths. The equations of motions 
are Xi = -{l/m)dH/dxi - irfLXi - FLit)) 5i^i - 
[irXn ^ Fii{t)) 5i_N ■ In these equations a, 6, Xeq, 
7 and m are constants and Fxit), K — L,R are 
Gaussian random forces that satisfy {FK{t)FK{0)) = 
'2lKkBTK5{t)/m. We take — 7(1 — x) and 7r = 7(1 + 
x) jIxI < Ij and study the ratio between AJ = J{Tl = 
Th;TR = T,) + J{Tl = T,;Tr = Th) and Jq = |J(x = 0)|, 
where Tc and Th denote low and high temperatures and 
where the heat current J is calculated as the average 
over sites, at steady state, of Ji = {—ii {dHi+i^i/dxi)) 
with = D (e-"(^'+i-^'-^='j) _ i)^ Fig. 1 shows 

that the intrinsic non-linearity of this model is enough 
to induce asymmetry in the thermal conduction of the 
asymmetrically coupled (7^ ^ ) bridge. 

Clearly, in addition to structural asymmetry, non- 
linear interactions are essential for rectifying behavior. 
In this paper we examine the rectifying properties of the 
simplest non-linear heat conductor: a two level system 
(TLS). The model investigated is a generalization of the 
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FIG. 1: The asymmetry in the thermal conduction plot- 
ted against x for ^ classical A'^-atom chain. The parame- 
ters used, D = 3.8/i/^ eV, a = IMv Xe9=1.538 A and 
m = 12/6.02 10~^^ g, are based on a standard model for the 
carbon-carbon force field in alkanes 16], for which u=l. Here 
we increase the system anharmonicity by taking v=fi. Full, 
dashed, dotted and dashed-dotted lines correspond to A'^^10, 
20, 40 and 80, respectively, with 7=50 ps-\ Tu = 300 K 
and Tc= K. (Inset) The temperature profile for the A''=80, 
X=0.5 case with Tl = T,; Tr = Th (fuU), Tl ^ Th ; Tr = 
(dashed). 



spin-boson model in which the TLS is coupled to two 
equilibrium boson baths maintained at different temper- 
atures. We study two variants of the model and show 
that if asymmetry is built into either one by employing, 
e.g. different spin-boson coupling strengths for the two 
baths, thermal rectification naturally sets in. 

The first variant of our spin-boson model is defined by 
the Hamiltonian 



H = Ec 



+ Ei\1){1\ + Hb + Hmb 



(3) 



Hb = Hl + Hr- Hk^Y. "^J^-W^ K^L,R (4) 

jeK 



H 



MB 



B|0)(l|+Bt|l)(0|; B = Bl+Br, (5) 



where a^, aj are boson creation and annihilation opera- 
tors associated with the phonon modes of the harmonic 
baths and Bk are bath operators. For a linear coupling 
model 



B 



K 



aj); K — L, R. 



(6) 

An important attribute of this model is the mu- 
tual independence of the transport processes at the two 
system-bath interfaces. It is convenient to regard the 
model (EJ-® as a special case of an N equally spaced 
states system with the nearest-neighbor coupling Hmb = 
^^JTj^^ -^71 (B|n — l)(n| -I- — In particular 



the N ^ oo limit corresponds to a harmonic oscilla- 
tor bridge connecting the baths. Eq. © corresponds 
in the latter case to the bilinear coupling model for 
the oscillator-baths interactions, Hmb = ^j^k ^j^j^^ 
where x is the coordinate of the bridge oscillator and 
aj = cij {21710)0)^ ^"^ , where m and ujq ^ Ei — Eq are the 
oscillator mass and frequency, respectively. 

The reduced dynamics of the TV-level system can be 
derived following standard procedures, e.g. the Redfield 
approximation |l8l | for the weak system-baths coupling 
limit. In the limit of fast dephasing the resulting kinetic 
equations for the state probabilities are 

Pn^ - {nkd + {n + l)fc„X„) P„ + nkuPn-i 

+ {n+l)kdXnPn+i; P-1^0, (7) 

where Xn ~ (5„.o ior the two level (?i=0,l) system and 
Xn=l for the harmonic oscillator (n — 0,..,oo ) case, 
and where kd = J^^dTe^'^"^ {B'< {t)B{0)) and ku = 
I-oo dTe~'"^°''{B{T)B'' (0)) . The average is over the baths 
thermal distributions, irrespective of the fact that it may 
involve two distributions of different temperatures |l9j. 
Specifying to the linear coupling model, and assuming no 
correlation between the thermal baths leads to the rates 



kn 



(8) 



with 



fci = rL(LJo)(l + nL{ujo)); kn = Tr{ujq){1 + urIluo)), 

(9) 



a^J{Lu-ujA; K^L,R. (10) 



2muj 



The heat conduction properties of this model are ob- 
tained from the steady state solution of Eqs. O with 
the rates given in Eqs. l(5|)- (|10() . For the harmonic model 
(N 00), putting Pn — 0, and searching a solution of 
the form P„ cx ?/" we get a quadratic equation for y whose 
physically acceptable solution is 



y 



(11) 



This leads to the normalized state populations P„ = 
?/"(l — y). The steady-state heat flux is obtained from 



J = Wo ^ ri {kjjPn - kj^Pn-ic 



(12) 



where positive sign indicates current going from left to 
right [23. Using Eqs. ® and JTTJ we find 



1 L + 1 i?. 



(13) 



This is a special case (with T(a;) = TiXRiJ^L + 
Tr)~^5{uj — Wo) consistent with our resonance energy 
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transfer assumption) |2j of Eq. Obviously no recti- 
fying behavior is obtained in this Umit. 

Next consider the two levels case, N=2. The two 
steady-state equations obtained from {Tj) yield 



Pi 



Re 



and the analog of Eq. H12() is 

J = uj^kR (Pi - Poe-''«"°) 
Using this with Eq. leads to 



(14) 
(15) 



J = Wq 



rL{l + 2nL)+TRil + 2nR)' 



which does have rectifying behavior. Indeed, defining 
the asymmetry parameter x such that = r(l — %) ; 
Tr = r(l + x) with -1 < X < 1 we find 



A J = Jin = n; Tr = T,) + Jin = r^; n = T;,) 
^ ^orx(l-x')K-nfl)^ 

(1 + TT-L + rt_R,)2 _ ;)(-2(^^ - ,^^)2 



(17) 



Eq. (inil implies that for smaU AT = Tl-Tr, |AJ| grows 
like AT^. Furthermore, noting that signiAJ) — signix) 
it follows from (|17|1 that the current is larger when the 
bridge links more strongly to colder reservoir than when 
it links more strongly to the hotter one. Figure 2 shows 
an example of this behavior. 




FIG. 2: Heat rectification by a TLS bridge in the lin- 
ear coupling model, Eqs. llHl- 11171 . The ratio A J/ Jo (with 
Jo = |J(x = 0)1 ) is plotted against the asymmetry param- 
eter X for several two-level spacings: ljo =0.025 eV (dashed); 
(.1^0=0.05 eV (full); ljo=0.075 eV (dotted). The temperatures 
are Th^400 K, Tc^SOO K. 

Next we consider another variant of the two-bath spin- 
boson model, taking the Hamiltonian to be 



H = Po|0)(0| + Pi|l)(l| + ■K),i|0)(l| + yi,o|l)(0| 
- ^ w,a]a,+ 2:,(ao,,|0)(0|+ai,,|l)(l|). (18) 

j£L,R j&L,R 



The L and R boson baths are again maintained at differ- 
ent temperatures n and T^. When Tl = Tr, Eq. lfTH|) 
represents a standard spin-boson Hamiltonian used, e.g, 
in the electron transfer problem. Using the small polaron 
transformation 22], H — UHU^^, leads to 



H = Eo 

jeL,R 



+ Pi|l)(l| + yo,i|0)(l|e''' + yi,o|i)(o|e-^" 

shift, (19) 



where U — UqUi 
A„., = (2c.3)-i/2 



Un 



(16) HsMft ^-il/2)j:^ 



expi-inn\n){n\), (n=0,l), 
(a]-aj) (if = L,i?), 

an.j and Q — Cli — CIq. The term 
4j0){0\+aljl){l\) maybe 



henceforth incorporated into the zero order energies. The 
Hamiltonian H19(l is similar to that defined in Eqs. ©- 
(O, except that the system-baths couplings appear as 
multiplicative rather than additive factors in the interac- 
tion term, implying non-separable transport at the two 
contacts. The dynamics is still readily handled. For small 
V (the " non-adiabatic limit" ) the Hamiltonian H19|l leads 
again to the rate equation Q with 

fcd= |K),il'C(c^o); K^\VoM''Ci~iUo), (20) 

where C(t^o) = JZc dte''^'>*Cit) and 



,i(0f(t)-0«(t))„-*(0f-^0«) 



L 

21) 



This may be evaluated explicitly to give 

Cit) = CLit)CRit); CKit) = expi-cl>Kit)) 



Aoj) [(l + 2n/f(wj)) 



(1 -I- riKiujj)) e~ 



(22) 



(23) 



Explicit expressions may be obtained using the short time 
approximation (valid for ^-^Oj)^ ^ 1 and/or 

at high temperature) whereupon (/)(t) is expanded in pow- 
ers of t keeping terms up to order t^. This leads to 



Ciu^) 



2-K 



-iujQ - E^j - EfiY 



where E^j 



(24) 

Xo,j)^Lu] i2nKiujj) + 1) 2kBTKE^j. Eqs. 

(|24|l provide an extension of the Marcus non-adiabatic 
rate expressions |2^ to the case of two reservoirs main- 
tained at different temperatures. E^j and E^j are the 
corresponding reorganization energies. 

Consider now the steady state heat current. The non- 
separability of the system-bath couplings makes the pro- 
cedure that leads to Eq. (|15|l unusable. Instead note that 
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Cl(wo) and Ci{(wo) are the rates affected by each ther- 
mal reservoir separately and that, from H22|l . C{u!o) = 
J^^dujCL{uJo—iLj)Cii{uj). The process |1) |0) in which 
the TLS looses energy luq can be therefore viewed as a 
combination of processes in which the system gives en- 
ergy uj (or gains it if < ) to the right bath and energy 
ujQ — uj to the left one, with probability Cl(clIo — ^)C'fl(a;). 
A similar analysis applies to the process |0) The 
heat flux calculated as the energy transferred per unit 
time into the right bath is therefore 



J = \Vo,l\' I dL0Lu[CR{L0)CL{i0Q - L0)Pl 

-Cr{^uo)Cl{-ujo+uj)PoI (25) 



where Pq — C {loq) / {C [loq 
Pq are the steady state 
tern is in state or 1, 
time approximation C{uj) 



-1/2 



exp 



-{u- 
Em{1 

which implies 
Ts =Tl + Tr. Using these 



to take Ej(^ 



- x); 
n2 - 



) + C{-uq)) and Fi = 1 - 
probabilities that the sys- 
respectively. In the short 
takes the form Ck{^) = 
/2Z?|j-]. It is convenient also 

E^j^Em{1+x) (IxI < 1), 
= 2kBEM{Ts - xAT) where 
relationships in l|25|l leads to 



J = 



2V^|Vo,il'(l 



AT 



3/2 I g2uJo/kB{Ts-x^Ty 



(26) 



Eq. H26|) again implies asymmetric heat conduction pro- 
vided symmetry is broken by taking x 0- This is shown 
in Fig. 13 where AJ/Jq is displayed against x- It is 
seen that the heat conduction asymmetry can be quite 
large, with its magnitude and sign strongly dependent 
on system parameters. When Em ^ wq the heat flux 
is dominated by the term e-('^o-2£^")V4feBBM (t,-xAT) 

that is bigger when AT is negative than when it is pos- 
itive, hence the negative asymmetry in AJ. The same 
behavior is seen in the opposite limit, Em ujq- How- 
ever, when 2Em ~ and wq ~ ksT, J is dominated 
by the term {ksEMiTs — x^T))^^^^ , implying positive 
asymmetry as seen in Fig. 13 

In Summary, while rectification of electronic current in 
molecular junctions is well known, heat flux rectification 
is a novel concept. Asymmetric anharmonic chains have 
this property. We have presented a simple heat rectify- 
ing model where anharmonicity stems from the dynamics 
of a two-levels system and asymmetry is introduced by 
different interaction strengths with the thermal baths. 
We have considered two cases that yield analytical solu- 
tions: a separable model with additive interactions with 
the bridge, and a non separable model. For both, the 
calculated heat current shows diode like behavior that 
depends on the junction characteristics. 

Asymmetric coupling to the two thermal reservoirs can 
be affected by different chemical bonding or by using 
phonon reservoirs with different Debye temperatures. Al- 
ternatively, in a more realistic molecular model, this may 




FIG. 3: Heat rectification of a TLS bridge with a non- 
separable coupling, Eqs. 1181 - 1261 . A J/ Jo is plotted against 
X for a;o=0.025 eV and for Em=0M2 eV (dashed), Em =0.38 
eV (full). The temperatures are 1^=400 K, Tc=300 K. 



result from asymmetric spatial organization of the molec- 
ular vibrational states. 

Heat rectification will be very useful in nanodevices, 
where efficient heat transfer away from the conductor 
center is crucial for proper functionality and stability. 
Similarly, directed energy flow in biomolecules such as 
proteins |24l | may play a role in controlling conforma- 
tional dynamics. From the theory perspective, the effi- 
ciency of rectification and its possible optimization are 
issues for future considerations. 

This research was supported by the Israel National Sci- 
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